Lattice simulations of QCD with \ib ^ versus phase quenched QCD 
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Previously published lattice results for QCD at /is 7^ are compared to analytic predictions for 
phase quenched QCD. We observe that the strength of the sign problem in QCD is linked directly 
to the position of the phase transition line for pion condensation in phase quenched QCD and that 
the number of terms needed in the Taylor expansion approach depends on the strength of the sign 
problem. Together this emphasizes the physical importance of the sign problem and helps to clarify 
the range over which the Taylor expansion approach is practically applicable. Finally, we observe 
that the positions of the endpoint of the first order chiral phase transition in the QCD phase diagram 
found in two successive computations by Fodor and Katz are both close to the position of the phase 
transition line for pion condensation in phase quenched QCD. 



INTRODUCTION 

Results from lattice simulations for unquenched QCD 
at non-zero baryon chemical potential give us valuable in- 
formation about the non-perturbative sector of strongly 
interacting matter. The central aim of the lattice sim- 
ulations at non-zero baryon chemical potential is to de- 
termine the phase diagram of unquenched QCD. The fo- 
cus at present is on the chiral transition as a function of 
temperature and rather small chemical potential. Pre- 
dictions for the slope of the chiral crossover line at zero 
chemical potential and for the position of the endpoint of 
the first order chiral phase transition have already been 
set forth. Lattice simulations at non-zero baryon chemi- 
cal potential are, however, different from lattice simula- 
tions at zero baryon chemical potential in that they have 
to deal with the sign problem. The term 'sign problem' is 
used to describe the numerical difficulties in doing Monte 
Carlo sampling on a non-positive weight. The presence of 
the sign problem is a direct consequence of the imbalance 
between quarks and anti-quarks one imposes in order to 
generate a non-zero baryon density. Moreover, the sign 
problem has several direct physical implications. One of 
the most dramatic can be seen if we compare to phase 
quenched QCD. 

Phase quenched QCD is different from QCD in that 
one takes the absolute value of the fcrmion determinant 
in the measure 

Z PQ = J [AA V ]\ detCA/y*, + M7o + m)|^e- 5 ™. (1) 

The phase diagram of this theory has an extended re- 
gion which is dominated by a Bose-Einstein condensate 
of pions. Returning to QCD by including the phase of 
the determinant wipes out this pion phase entirely. In 
order to deal with the phase of the determinant in lattice 
simulations of QCD three main approaches have been 
pursued: 1) the Taylor expansion approach 0, 0, & @ > 
2) the reweighting approach 0, 0, W, SJl^nd 3) the 
imaginary chemical potential approach [lfj, |llj, |l2| . 
In this paper we address some of the issues involved 



in the interpretation of existing lattice data for QCD at 
non-zero baryon chemical potential obtained with these 
approaches. We will do so by comparing previously pub- 
lished lattice data for QCD at non-zero baryon chemical 
potential to phase quenched QCD. First of all we replot a 
lattice measurement from p] of the variance of the phase 
of the determinant in units where we can compare to an- 
alytic predictions for phase quenched QCD. We observe 
that the contour lines for the variance of the phase of the 
determinant are aligned with the phase transition line 
for pion condensation in phase quenched QCD. More- 
over, the distance between the contour lines decreases as 
the values of n and T approach the region where phase 
quenched QCD is in the pion phase. Next, we wish to 
verify, as stated in Q , that the order of the Taylor expan- 
sion needed depends on the strength of the sign problem. 
In order to do so we consider the 6th order Taylor ex- 
pansion of the quark number susceptibility published in 
[5j- In agreement with the statement of [J] we find that 
the 6th order term in the expansion becomes important 
when the variance of the phase exceeds a certain value. 
Together with the first observation we conclude that it 
will be exceedingly hard to make predictions for QCD 
at non-zero baryon chemical potential using the Taylor 
expansion approach when the values of /i and T are such 
that pion condensation occurs in phase quenched QCD. 

Finally, we consider the values of [i and T at the end- 
point of the first order chiral phase transition as found 
using the reweighting method in || and || . After rescal- 
ing the coordinates for the endpoint we observe that both 
predictions occur at values of fj, and T which are close 
to the phase boundary for pion condensation in phase 
quenched QCD. 



PION CONDENSATION IN PHASE QUENCHED 

QCD 

At zero temperature the critical chemical potential for 
pion condensation in phase quenched QCD is fi — 
and at small temperature the critical chemical potential 
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can be evaluated in chiral perturbation theory [13|. The 
critical temperature found at 1-loop order corresponds 
to the semi-classical result for the critical temperature of 
Bose-Einstein condensation in a dilute massive Bose gas. 
The general determination of this phase transition can be 
evaluated directly by means of Monte Carlo simulations 
of phase quenched QCD and such studies are currently 
taking place ^ij. One can also estimate the position 
of the phase transition using phenomenological models. 
Probably the simplest prediction for the critical temper- 
ature in phase quenched QCD is obtained from the ran- 
dom matrix model considered in [lij . This mean field 
ap proa ch leads to a critical temperature (see eq (5.48) 
of [15J) which depends both on the chemical potential, 
the quark mass and the phenomenological constant in- 
troduced in the model. However, by taking the quark 
mass to zero while keeping the ratio of the pion mass 
to the chemical potential fixed the dimensionful constant 
drops out and the result for the critical temperature for 
pion condensation in phase quenched QCD is (Tq is the 
pseudo critical temperature for chiral symmetry breaking 
at /z = 0) 

In the range < fi < this simple result a gree s 

well with the predictions from the N JL model [lJl Qj| , 
the prediction from strong coupling QCD 0], as well 
as the predictions from a random matrix model with all 
Matsubara frequencies included [IT]]. At larger values of 
the chemical potential the critical temperature obtained 
in these models drops down to zero again as a result of 
saturation. 

The fi dependence of the critical temperature given in 
is also consistent with the picture e merg ing from the 
lattice studies of phase quenched QCD 1141. These lat- 
tice simulations and analytical arguments |l3l llif how- 
ever suggest that the order of the pion phase transi- 
tion changes from 2nd order to 1st order with increasing 
chemical potential. This aspect is not reproduced by the 
random matrix model. 

THE VARIANCE OF THE PHASE 

A direct measure of the sign problem in QCD is the 
variance of the phase, 9, of the fermion determinant. 
The phase was measured in 0, l2p| and a contour di- 
agram with lines of constant variance, \J (9 2 ) — (9) 2 , 
in the (T/Tq, /i/T)-plane was given. Below we re- 
plot the data from as a contour diagram in the 
(2/i/m T , T/Tq) -plane, To again being the pseudo-critical 
temperature for the crossover at /i = 0. The value of 
the pion mass used in order to convert the plot from 
5] is m^/To — 3.58. Shown are contour lines for 




0.5 1 1.5 2 



2\x/m % 

FIG. 1: In the (2^/m 7r ,r/T )-plane we plot 

• In red: Lines of constant variance 

y/W)~W horn g. 
From the upper left the values increase from 7r/4 to 2ir in 
steps of 7r/4. Lines of higher variance are not available. 

• In black: The critical temperature as a function of the chem- 
ical potential for pion condensation in phase quenched QCD 
as estimated in @. 

V( 62 ) - W = t/ 4 » 3tt/4, . . . , 2vr, the value increas- 
ing with the chemical potential. Note that contour lines 
for y/ (9 2 ) — (9) 2 > 2ir are not displayed. In addition to 
the contour lines the prediction (j2J for the critical tem- 
perature for pion condensation in phase quenched QCD 
is shown. The contour lines are seen to be aligned with 
the phase transition line and the distance between the 
contour lines decreases as the phase transition line is ap- 
proached. This observation is consistent with the fact 
that the presence of the phase of the determinant wipes 
out the pion phase transition. In particular, it is natural 
to expect that ^ {9 2 ) - (9) 2 » 1 in the region where 
phase quenched QCD enters the pion phase, even with 
a moderate volume. At zero temperature this has been 
verified analytically [2ll : the eigenvalue density of the 
Dirac operator in QCD with /i ^ becomes highly os- 
cillating for /i > m T /2 and it is these oscillations which 
insure that the pion phase is avoided in QCD 23]. Note 
that the contour lines shown in figure ^ will shift toward 
lower values of n and higher values of T with increasing 
volume (for a discussion of the volume dependence of the 
variance see |2oj]h 



TAYLOR EXPANSIONS 

One of the approaches to QCD at non-zero quark num- 
ber chemical potential which are being pursued 
is based on a Taylor expansion in \ijT around // = 0. 
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FIG. 2: Figure from y showing the quark number suscep- 
tibility for fixed values of fi/T as a function of temperature. 
Full lines give the result of the 6th order Taylor expansion 
while the dashed lines give the results at 4th order. 

First, the forms of the traces contributing to an operator 
at a given order in the chemical potential are determined 
analytically and then these traces are evaluated by stan- 
dard Monte Carlo methods at /i = 0. Traces contributing 
to order 2n in the Taylor expansion of an extensive ob- 
servable are typically of order V n and these terms must 
combine to a result linear in the volume 4]. This makes 
it difficult to control numerical errors in a high order ex- 
pansion. 

In p| it is stated that the order of terms needed in 
the Taylor expansion depends on the severity of the sign 
problem. Here we verify this statement using the quark 
number susceptibility determined in || to 6th order in 
li. The relevant figure from || is given in figure El for 
convenience. The quark number susceptibility is shown 
as a function of temperature for fixed values of li/T as 
obtained from the 6th order Taylor expansion. For the 
three largest values of li/T the result from the 4th order 
Taylor expansion are indicated by the dashed lines. From 
this plot we estimate the values of T/Tq below which the 
6th order term is important. Using simply the magnitude 
of the deviation in the units of figure for /u/T = 1.0 we 
find T/Tq ~ 1.20, for fx/T = 0.8 we find T/Tq ~ 1.10, 
and for ll/T = 0.6 we find T/Tq ~ 1.05. In figure we 
plot the contour lines for the variance of 9 from |5j in the 
(2/x/to w , T/Tq) plane, this time together with the 3 lines 
determined by li/T = 1.0, 0.8, and 0.6. On the lines of 
constant li/T we have just estimated the values of T/Tq 
below which the 6th order term in the Taylor expansion 
of the quark number susceptibility becomes important. 
These values are indicated by the three crosses in figure 
|5J We observe that the 6th order term in the determi- 
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FIG. 3: In the (2^/m^,T/To)-plane we plot 

• Red lines: Lines of constant yj {8 2 } — (8) 2 from From 
the upper left the values increase from 7r/4 to 2n in steps of 
tt/4. 

• Black lines: Lines with fixed ratio of (i/T, the values shown 
are fi/T = 0.6, 0.8, and 1.0. 

• Black crosses: Marks the point on the black lines below 
which the 6th order term in the Taylor expansion of the quark 
number susceptibility in [f| is important. 

nation of the quark number susceptibility becomes im- 
portant for values of T and li where the variance of the 
phase is larger than 7r/4. That is, the 4th order Tay- 
lor expansion of the quark number susceptibility in |5j is 
only appropriate to the left of the leftmost contour line 
in [5] Given the sufficient numerical accuracy to handle 
the delicate cancellations of the terms of order V 6 , the 
6th order Taylor expansion of the quark number suscep- 
tibility approach may be adequate somewhat to the right 
of the y/ (8 2 ) - (0) 2 = 7r/4 contour line. 

Combined with the observations made from figure^wc 
conclude that even for small volumes it will be extremely 
difficult to make predictions for observables in QCD using 
the Taylor expansion approach when the values of ll and 
T are such that phase quenched QCD is in the pion phase. 

Certainly, our estimate above only serves as a first nu- 
merical test of the proposition that the number of terms 
needed in the Taylor expansion depends on the strength 
of the sign problem. A dedicated series of measurements 
of several observables on the lattice are needed in order 
to settle this firmly. 

CRITICAL ENDPOINT 

In this section we compare the positions of the end- 
points of the first order chiral phase transition as deter- 
mined in Q and to the estimated position of the pion 
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and the eigenvalue density is 

Ppq(z, to, T, h) = I y S 2 { 



(5) 
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FIG. 4: In the (2^/m 7r ,T/T )-plane we plot 

• Red squares: The position of the critical endpoint as com- 
puted in (lower point) and in [IJ (upper point). 

• Black line: The critical temperature as a function of the 
chemical potential for pion condensation in phase quenched 
QCD as estimated in @. 



phase transition in the phase quenched theory 1 . In order 
to facilitate this comparison we again consider the plane 
(2/z/to^, T/Tq) and plot in figure 0] simultaneously the 
critical endpoints determined for QCD in 0, and the 
position of the pion phase transition for phase quenched 
QCD given in (J2J. The values of the pion mass used in 
order to make the rescaling was r?v = 294MeV for 
and = 145MeV for 0. 

We observe that the endpoints found in [8j and |9j 
both are remarkably close to the phase transition line 
estimated in for phase quenched QCD. (The fact that 
the chemical potential at the endpoint scales with m n was 
noted in 0].) 

The phase transition into the pion phase occurring in 
phase quenched QCD has a very simple manifestation in 
the eigenvalue spectrum of the Dirac operator: it occurs 
when the density of eigenvalues of the Dirac operator 
reaches the quark mass (see e.g. [24| or H2,[2^]). That 
is, the critical chemical potential and temperature are 
determined by 

Ppq(z = to, to, T c ,p c ) 0, (3) 

where the eigenvalues are given by 

(A;77J + M7o)^n = z n^n (4) 



PQ 



1 I thank Misha Stephanov for suggesting this comparison. 



For a lattice simulation involving the full determinant of 
the QCD Dirac operator this leads to substantial prob- 
lems: For /x larger than the critical /x for pion conden- 
sation in phase quenched QCD the eigenvalues z n of 
Dnii-i + M7o lie arbitrarily close to —to. The complex 
phase of (z„ + m) therefore is extremely sensitive to the 
gauge field configuration and this manifests itself in the 
fluctuation of the phase of the fermion determinant. For 
this reason the reweighting methods are exceedingly del- 
icate for values of [i and T where phase quenched QCD 
is in the pion phase. In the light of this, it is unfortunate 
that the predictions for the critical endpoint in both 
and are so close to the predicted phase transition for 
phase quenched QCD. 



CONCLUSIONS 

In order to emphasize the physical importance of the 
sign problem for present lattice calculations at non-zero 
baryon chemical potential we have replotted previously 
published lattice data in units where a direct comparison 
to phase quenched QCD is possible. We have observed 
that the seriousness of the sign problem, as measured by 
the variance of the phase of the fermion determinant, de- 
pends on the distance in the (2ii/m T , T/T )-plane from 
pion phase transition in phase quenched QCD. The esti- 
mate used for the position of the phase transition in phase 
quenched QCD is independent of any phenomenological 
parameters. It is obtained as the limit of small quark 
masses with [i / 'TO^-fixed of the result from a random ma- 
trix model. While this result is in reasonable agreement 
with other phenomenological models one should only take 
it as a rough guideline for the position of the pion phase 
transition in phase quenched QCD. Lattice simulations 
of phase quenched QCD can measure the position as well 
as the order of the phase transition and thus allow for 
a direct comparison between QCD and phase quenched 
QCD. 

Using the result for the Taylor expansion of the quark 
number susceptibility in |5| we have verified that the 
number of terms needed in the Taylor expansion ap- 
proach depends on the seriousness of the sign problem 
as measured by the variance of the phase of the fermion 
determinant. The need for a Taylor expansion to 2rdh 
order post a demand for sufficient numerical accuracy to 
realize cancellations between terms of order V n in or- 
der to obtain an expectation value of order V @]. Thus 
the sign problem sets the practical upper bound on the 
applicability of the Taylor expansion approach. In par- 
ticular it will be extremely computationally demanding 
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to determine observables in QCD by means of a Taylor 
expansion if the values of fj, and T are such that phase 
quenched QCD is in the pion phase. 

The need for a high order Taylor expansion to pene- 
trate regions of the (/i,T)-plane where the sign problem is 
strong also sets a practical limit on the imaginary chem- 
ical potential approach. The analytic continuation from 
imaginary to real chemical potential is carried out by 
fitting a polynomial in fi and the order of the polyno- 
mial needed corresponds to that of the Taylor expansion. 
We have checked that no data points have been reported 
in the literature with this method for values of \x and 
T where phase quenched QCD is in the pion condensed 
phase. 

Finally, we have observed that the location of the end- 
point of the first order chiral phase transition determined 
using the reweighting method in [8| as well as in |9| falls 
surprisingly close to the location of the pion phase transi- 
tion in the phase quenched theory. Despite attempts, no 
analytical argument to date suggests that this should be 
the case. For values of T and [i where the phase quenched 
theory enters the pion phase the reweighting approach 
is expected to be extremely delicate. Hence, one may 
worry that the critical point in || as well as in || are 
manifestations of the numerical difficulties encountered 
rather than a true physical effects. Certainly it would be 
interesting if by an appropriate choice of quark masses 
one could separate the prediction for the critical point in 
QCD from the phase boundary in phase quenched QCD. 
As suggested in ^l) fine tuning the quark masses close 
to critical value for which the chiral transition at fi = 
becomes first order may cause the endpoint to move to 
smaller values of 2u./m w . 
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